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addenda 


1.  On  9  June  1982 ,  the  Defense  Mapping  Agency  issued  upgraded 
WGS-72  coordinates  for  the  Loran-C  stations.  This  update  includes 
positions  measured  to  the  nearest  0.001",  the  *7930  Northwest 
Pacific  Reconf iguation  and  the  addition  of  the  5970  Commando  Lion 
Chain.  The  upgraded  coordinates  are  listed  on  the  following  page. 

2.  The  HP-41CV  program  may  be  modified  to  display  the  station 
coordinates  to  the  nearest  0.001"  by  changing  line  953  to  F1X2. 
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ABSTRACT 


This  report  presents  an  algorithm  and  HP-41CV  programs  for 
position  determination  with  Loran-C  chains.  Additional 
computational  routines  include  the  ability  to  calibrate  Loran 
station  triplet  data  to  a  known  benchmark  and  ITD’s  (Indicated 
Time  Delay's),  predict  ITD's  at  given  positions,  compute  the 
geodesic  (similiar  to  great  circle)  bearing  and  distance  from  a 
fix  to  the  destination  and  to  compute  the  geodesic  bearing  and 
distance  from  any  one  location  to  another.  Utility  routines  allow 
the  user  to  transfer  station  pair  data  between  the  HP-41CV  and 
magnetic  cards,  magnetic  tape  and  an  extended  function/  memory 
module. 


A.  Introduction 

The  Loran  system  is  a  radio  aid  to  navigation  which  utilizes 
the  principle  of  hyperbolic  fixing.  The  locus  of  points  for  which 
the  difference  in  arrival  time  of  synchronized  signals  from  a  pair 
of  transmitters  is  constant  determines  a  hyperbolic  line  of 
positions.  The  intersection  of  two  hyperbolic  lines  of  position 
from  two  pairs  of  stations  determines  position  or  a  hyperbolic 
fix.  That  two  pairs  of  stations  are  required  for  a  fix  does  not 
necessarily  mean  that  there  are  four  separate  stations,  for  one 
station  of  one  pair  may  be  colocated  with  one  station  of  the  other 
pair  forming  a  Loran  triplet  (Figure  1).  Triplets  may  be  joined 
"end-to-end"  by  station  colocation  to  form  a  Loran  chain  (Figure 
2).  Loran  chains  are  common  on  both  the  East  and  West  coasts  of 
the  North  American  continent. 

The  present  day  Loran-C  operates  at  100-kHz  and  is  in  use  in 
the  Atlantic,  Pacific  and  Mediterranean  areas.  The  computational 
algorithm  and  programs  described  herein  can  be  used  for  position 
determination  with  Loran-C  triplets.  Further  information  on  the 
history,  development  and  operation  of  the  Loran  systems  may  be 
found  in  References  1  and  2. 
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(a)  Colocated  Master  Stations 


(b)  Colocated  Slave  Stations 


(c)  Colocated  Master  and  Slave 

Figure  1.  Loran  Triplets 


B.  user  Routines 


1.  Computational  Routines: 

FI  -  The  Fixing  routine  is  the  main  program  for  calculating 
a  Loran-C  fix  from  indicated  time  delays. 

AS  -  The  Alternate  Solution  routine  will  allow  the  second 

Loran  fix  solution  to  be  computed.  This  routine  toggles 
Flag  3  so  that  on  subsequent  fixes  the  FI  routine  will 
calculate  the  alternate  solution. 

DN  -  The  destination  routine  stores  the  latitude  and 
longitude  of  a  fixed  destination. 

HD  -  Computes  the  Heading  and  distance  from  the  current  fix 
to  the  destination  stored  by  the  DN  routine. 

2.  Manual  Mode  Routines: 

MI  -  Manual  Input  allows  station  data  to  be  input  and  stored 
via  the  calculator  keyboard. 

fid  -  ficho  data  is  a  utility  routine  for  validating  station 
triplet  information  stored  in  the  calculator. 

3.  Card  Reader  Routines: 

CS  -  Card  Store  records  station  data  onto  magnetic  cards. 

CR  -  Sard  Read  inputs  station  data  from  magnetic  cards. 

CE  -  Card  £cho  is  a  utility  routine  for  validating  station 
information  stored  on  data  cards. 

4.  Extended  Memory  Routines: 

AS  ~  XMEM  Store  records  station  data  onto  the  extended 
memory  module. 

XR  -  XMEM  Recall  inputs  station  data  from  the  extended 
memory  module. 
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HD  -  XMEM  Delete  erases  station  data  from  the  extended 


memory  module. 

5.  Tape  Cassette  Routines: 


is  - 

Tape 

Store 

records 

station 

data 

onto 

the 

tape 

cassette 

Tfi  - 

Jape 

fiecall 

inputs 

station 

data 

from 

the 

tape 

cassette 

ID  - 

Tape 

Delete 

erases 

station 

data 

from 

the 

tape 

cassette 

6.  Utility  Routines: 

D£L  -  This  routine  is  similiar  to  the  HD  option  except  that 

it  computes  the  heading  and  distance  from  any  origin  to 
any  destination. 

-  PR  is  used  to  PRedict  the  station  ITD's  that  will  be 
received  at  a  given  latitude  and  longitude. 

DA  -  CA  is  the  Calibration  option.  Given  the  latitude  and 
longitude  of  a  known  position  and  the  indicated  time 
delays  from  a  Loran-C  triplet,  the  stored  station  data 
are  modified  so  that  the  FI  routine  (or  AS)  will 
compute  the  known  position  from  the  same  time  delays. 

SH  -  SHitch  data  swaps  the  data  of  the  two  Loran  stations 
stored  in  memory. 

Note:  There  are  no  specific  routines  that  relate  to  HP-41CV 
printer  operations.  However,  all  input  and  output  will  be 
recorded  on  a  printer  if  one  is  attached. 
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c.  Recommended  Eey  Assignments 

It  is  recommended  that  the  following  HP-41CV  user  key 
assignments  be  made  and  recorded  onto  the  program  cards  if  the 
program  is  first  prepared  manually: 


CA  ED 


1 

CS 

CR 

CE 

TD 

MI 

sw 

TS 

TR 

DN 

XD 

PR 

DH 

XS 

XR 

HD 

AS 

FI 

u 

The  El  function  is  placed  on  the  ENTER  key  as  a  reminder  that 
ENTER  should  npt  be  used  for  data  entry.  To  make  the  user  key 
assignments  on  the  HP-41CV,  refer  to  the  ASN  function  in  the 
HP-41 CV  Owner's  Handbook  and  Programming  Guide. 
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D.  User  Instructions  with  Examples 
1.  Manual  Input  Routine  -  MI 

This  routine  can  be  used  to  enter  station  data  (Appendix 
C)  to  prepare  data  to  be  transferred  to  cards,  extended  memory, 
or  tape  cassette.  It  can  also  be  used  to  enter  station  data 
manually  if  alternate  data  storage  media  are  not  available. 

Example  using  9940W  and  9940Y.  (Note:  The  notation  <9>, 
for  example,  means  that  you  must  press  the  gold  shift  key  and 
then  the  9  key.) 


See 

Key,  in 

Press 

1. 

XEQ"MI" 

2.  ID? 

<9><9><4><0>W 

R/S 

3.  CODE  DELAY? 

11000 

R/S 

4.  MS  LAT? 

39.330662 

R/S 

5.  MS  LON? 

118.495637 

R/S 

6.  SS  LAT? 

47.034799 

R/S 

7.  SS  LON? 

119.443953 

R/S 

If  desired,  these 

station  parameters  can  now  be 

transferred  to  card,  extended  memory  or 

tape  using  the  CS,  XS, 

or  TS  routines,  respectively.  Otherwise 

,  repeat  the  steps  above 

with  the  9940Y  station 

data.  These  two 

stations  will  be  used  in 

the  remaining  examples. 

See 

in 

PLfiSfi 

1. 

XEQ"MI " 

2.  ID? 

<9><9><4><0>Y 

R/S 

3.  CODE  DELAY? 

40000 

R/S 

4.  MS  LAT? 

39.330662 

R/S 
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5.  MS  LON? 


118.495637 


R/S 

6.  SS  LAT?  35.191818  R/S 

7.  SS  LON?  114.481743  R/S 

8.  NEXT  OPTION? 

At  this  point  the  station  data  for  9940W  and  9940Y  are 
stored  in  the  calculator.  The  9940Y  data  can  be  transferred  to 
card,  extended  memory  or  tape  cassette  using  the  CS,  XS  or  TS 
routines,  respectively. 

[Advanced  User  Note:  The  first  action  taken  in  the  MI 
routine  is  to  XEQ"SW".  SW  is  the  data  swap  routine  which 
exchanges  to  content  of  R16  -  R25  with  R26  -  R35 .  The  incoming 
data  are  then  stored  in  the  R16  -  R25  registers.  The  CS,  XS  and 
TS  routines  transfer  the  content  of  R16  -  R25  to  card,  extended 
memory  or  tape  cassette,  respectively.  If  desired,  the  content 
of  these  registers  can  be  swapped  once  more  by  using  the  SW 
utility  routine.] 
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.  Echo  Data  Routine  -  ED 

This  routine  allows  the  user  to  review  the  station  data 
resident  in  the  calculator. 

Example:  Load  the  station  pairs  9940W  and  9940Y  using 
either  the  MI,  CR,  XR  or  TR  routines. 


See 

Key  in 

Press 

1. 

XEQ"ED" 

2.  ID:  9940W 

R/S 

3.  CD:  11000 

R/S 

4.  MLT:  39.330662 

R/S 

5.  MLN:  118.495637 

R/S 

6.  SLT:  47.034799 

R/S 

7.  SLN:  119.443953 

R/S 

8.  BL:  2796.903 

R/S 

9.  ID:  9940 Y 

R/S 

10.  CD:  40000 

R/S 

11.  MLT:  39.330662 

R/S 

12.  MLN:  118.495637 

R/S 

13.  SLT:  35.191818 

R/S 

14.  SLN:  114.481743 

R/S 

15.  BL:  1967.302 

R/S 

16.  NEXT  OPTION? 

Notation:  CD  =  coding  delay,  M  *  master,  S  =  slave,  LT  = 
latitude,  LN  =  longitude  and  BL  is  the  station  pair  baseline 
plus  the  secondary  phase  correction  in  microseconds. 
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3a.  Card  Store  Routine  -  CS 

With  the  card  reader  attached,  station  data  (in  R16  - 
R25),  which  has  been  input  using  the  MI,  XR  or  TR  routine,  can 
be  transferred  to  magnetic  card  using  the  XEQ"CS"  command  or  by 
pressing  the  appropriate  user  defined  key. 

Example  using  9940W. 

Key  in  £X£££ 

1.  XEQ"PD " 

2.  WRITE:  9940W  (Pass  a  blank  card  through  the  card 

reader.  Label  the  card  track  "9940W") 

3.  NEXT  OPTION? 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  prepare  a  card  for  the  9940Y  station  pair. 


3b.  Card  Read  Routine  -  CR 

With  the  card  reader  attached,  XEQ"CR"  or  press  the 
appropriate  user  defined  key.  This  routine  can  be  used  to  input 
the  data  for  two  station  pairs,  which  must  form  a  triplet. 
Example  using  9940W  and  9940Y. 

See  Key  in  Press 


1.  XEQ"CR" 

2.  1ST  CARD  (Pass  the  data  card  for  9940W  through 

the  card  reader.) 

3.  2ND  CARD  (Pass  the  data  card  for  99440Y  through 

the  card  reader.) 

4.  NEXT  OPTION? 
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3c.  Card  Echo  .Routing--  CE 


With  the  card  reader  attached,  XEQ"CE"  or  press  the 
appropriate  user  defined  key.  This  routine  is  used  to  validate 
the  content  of  data  cards  against  the  table  in  Appendix  C. 
Example  using  9940W. 

Key  in  Press 

1.  XEQ"CE" 


2.  STA.  CARD  (Pass  one  side  of  a  data  card  [9940W] 
through  the  card  reader) . 


3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 


ID:  9940W 

R/S 

CD:  11000 

R/S 

MLT:  39.330662 

R/S 

MLN:  118.495637 

R/S 

SLT:  47.034799 

R/S 

SLN:  119.443953 

R/S 

BL:  2796.903 

R/S 

NEXT  OPTION? 

Notation:  CD  =  coding  delay,  M  =  master,  S  *  slave,  LT  = 
latitude,  LN  =  longitude  and  BL  is  the  station  pair  baseline 
plus  the  secondary  phase  correction  in  microseconds. 
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4a.  Store  Data  in  Extended  Memory  -  XS 

With  the  extended  memory  module  in  the  HP-41CV,  station 
data  (in  R16  -  R25),  which  has  been  input  using  the  MI,  CR,  XR 
or  TR  routine,  can  be  transferred  to  the  module  using  the 
XEQ"XS"  command  or  by  pressing  the  appropriate  user  defined 
key.  Example  using  9940W. 

See  Key  in  Press 

1.  XEQ"XS" 

2.  NEXT  OPTION? 

Should  the  station  pair  already  be  in  extended  memory,  the 
message  DUP  FL  (duplicate  file)  will  be  displayed.  If  needed, 
the  duplicate  file  may  be  erased  using  the  XD  routine. 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  store  the  9940Y  station  pair. 

4b.  Recall  Data,  f ism  JE&tended  Memory  -  XR 

With  the  extended  memory  module  installed,  XEQ"XR"  or 
press  the  appropriate  user  defined  key  to  input  the  data  for  a 
station  pair. 

Example  using  9940W.  (Note:  The  notation  <9>,  for  example, 
means  that  you  must  press  the  gold  shift  key  and  then  the  9 
key.) 

See  Key  in  Press 

1.  XEQ"XR" 

2.  ID?  <9><9><4><0>W  R/S 

3.  NEXT  OPTION? 
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4c 


XEQ"XD"  or  press  the  appropriate  user  defined  key  to 
delete  the  specific  station  pair  data  from  the  extended  memory 
module. 

Example  using  9940W. 

See  Key  in  Press 

1.  XEQ"XD" 

2.  ID?  <9><9><4><0>W  R/S 

3.  NEXT  OPTION? 

The  message  FL  NOT  FOUND  will  be  displayed  if  the  file  you 
wish  to  delete  in  not  in  the  extended  memory. 

Note;  The  extended  functions/  memory  module  will 
accommodate  the  data  for  11  station  pairs.  The  extended  memory 
module  will  accommodate  the  data  for  an  additional  22  station 
pairs. 
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5a.  stoie.  Data  in  .the.  Tape- Caaselte.  -  is 


With  the  tape  cassette  attached  to  the  HP-41CV,  station 
data  (in  R16  -  R25),  which  has  been  input  using  the  MI,  CR,  XR 
or  TR  routine,  can  be  transferred  to  the  tape  using  the  XEQ*TS" 
command  or  by  pressing  the  appropriate  user  defined  key. 

Example  using  9940W. 

£££  Key  -iQ  Eneas 

1.  XEQ"TS" 

2.  NEXT  OPTION? 

Should  the  station  pair  already  be  in  extended  memory,  the 
message  DUP  FL  NAME  (duplicate  file  name)  will  be  displayed.  If 
needed,  the  duplicate  file  may  be  erased  using  the  TD  routine. 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  store  the  9940Y  station  pair. 

5b.  Recall.  Data.  fiom.  .the  .lapfi.  Caaae.tte  -TR 

With  the  tape  cassette  attached,  XEQ"TR"  or  press  the 
appropriate  user  defined  key  to  input  the  data  for  a  station 
pair . 

Example  using  9940W.  (Note:  The  notation  <9>,  for  example, 
means  that  you  must  press  the  gold  shift  key  and  then  the  9 
key.) 


&&£ 

Kfii. JLn 

Etess 

1. 

XEQ"TR" 

2. 

ID? 

<9><9><4><0>W 

R/S 

3. 

NEXT  OPTION? 
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5c.  Delete  Data  from  the  Tape  Cassette  -  TD 

XEQ"TD"  or  press  the  appropriate  user  defined  key  to 
delete  the  specific  station  pair  data  from  the  tape  cassette. 
Example  using  9940W. 

fi££  Key  in  Press 

1.  XEQ"TD" 

2.  ID?  <9><9><4><0>W  R/S 

3.  NEXT  OPTION? 

The  message  FL  NOT  FOUND  will  be  displayed  if  the  file  to  be 
deleted  is  not  on  the  tape. 
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6.  Loran-C  Fixing  Routines  FI  and  AS 


Given  the  indicated  time  delay  (ITD)  from  two  station 
pairs  which  form  a  triplet,  a  Loran-C  fix  is  obtained. 


Example : 

Load  9940W 

and  9940Y 

into  the  calculator 

the 

MI, 

CR,  XR  or  TR  routine.  The 

ITD  on  9940W  is  16019 

microseconds 

and  the  ITD 

of  9940Y 

is  42585  microseconds 

are 

you? 

See 

Key.  Id 

1. 

XEQ"FI " 

2. 

ITD: 

9940W 

16019 

R/S 

3. 

ITD: 

9940  Y 

42585 

R/S 

4. 

LAT: 

39.1419 

R/S 

5. 

LON: 

115.5052 

R/S 

6. 

NEXT 

OPTION? 

Since  you  are  on  a  boat,  you  know  that  you  cannot  be  in 
central  Nevada  at  39dl4'19"  North  and  115d50’52"  West.  Every 
Loran-C  fix  has  two  solutions,  so  in  this  case  you  must  use  the 
alternate  solution. 


Sae  Key. .in 

7.  XEQ"AS" 

8.  LAT:  35.0001 

9.  LON:  125.0009 

10.  NEXT  OPTION? 


EX£fi£ 

R/S 

R/S 


This  is  the  proper  solution  at  almost  exactly  35  degrees 
North  and  125  degrees  West.  Note  that  annunciator  3  (Flag  3) 
shows  in  the  display  indicating  the  alternate  solution.  If  you 
should  now  repeat  from  Step  1,  you  will  obtain  the  proper 
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solution  inuned lately. 

The  message  "E:  NO  TRIPLET"  will  appear  following  Step  1 
if  the  data  do  not  comprise  a  valid  triplet.  The  latitudes  and 
longitudes  of  each  station  pair  at  the  vertex  must  agree 
exactly.  Should  this  error  occur,  use  the  ED  routine  to  review 
the  resident  station  data. 

The  message  "E:  ITD  ERROR"  will  appear  following  Step  2  or 
3  indicating  that  the  ITD  you  keyed  in  is  inconsistent  with  the 
station  parameters.  Press  R/S  to  be  requeried  for  the  ITD. 
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7a 


If  you  know  the  latitude  and  longitude  of  your 
destination,  you  may  key  these  in  and  then  see  how  far  your  fix 
is  from  your  destination  and  what  the  geodesic  heading 
(similiar  to  great  circle  heading)  is  to  your  destination. 

Example:  Your  destination  is  Moss  Landing  at  about  36d48'N 
and  121d47'W.  Your  current  fix  is  35dN  and  125dW  (see  the  FI-SA 
example) .  First,  key  in  your  destination. 


See 

Key  in 

Press 

1. 

XEQBDN" 

2.  DEST  LAT? 

36.48 

R/S 

3.  DEST  LON? 

121.47 

R/S 

4.  NEXT  OPTION? 

R/S 

The  destination  is  now  stored  in  the  calculator  and  will  remain 
unchanged  until  you  use  either  the  DN  or  DH  options.  Also,  the 
latest  fix  is  stored  and  will  remain  unchanged  until  you  use 
either  the  FI,  AS  or  DH  options.  Now,  find  the  distance  and 
bearing  from  the  latest  fix  (see  the  FI-AS  example)  to  Moss 
Landing. 

£££  Key  in  Press 

1.  XEQ"HD" 

2.  N.MI:  190.38  R/S 

3.  BRG:  54.3411  R/S 

4.  NEXT  OPTION?  R/S 

So  the  distance  to  Moss  Landing  is  190.38  nautical  miles 
at  a  heading  of  54d34'll". 
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7b 


Given  the  latitude  and  longitude  of  an  origin  and 
destination,  this  routine  will  find  the  distance  and  heading 
from  one  to  the  other. 

Example:  How  far,  and  in  what  direction,  is  Corvallis, 
Oregon  (44d34'N,  123dl6'W)  from  Cupertino,  California  (37dl9N, 
122d02 *W) ? 


See 

Key  Id 

XEQ"DH" 

ORIG  LAT? 

37.19 

R/S 

ORIG  LON? 

122.02 

R/S 

DEST  LAT? 

44.34 

R/S 

DEST  LON? 

123.16 

R/S 

N.MI :  438.32 

R/S 

BRG:  353.0259 

R/S 

NEXT  OPTION? 

Thus  the  distance  is  438.32  nautical  miles  and  the 
direction  of  Corvallis  from  Cupertino  is  353d02'59". 
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8.  ITD  Prediction  Routine  -  PR 

As  an  aid  to  identification,  this  routine  will  allow  the 
user  to  determine  what  ITD's  should  be  received  at  a  given 
location. 

Example:  Suppose  that  you  know  you  are  somewhere  near 
latitude  35  North  and  longitude  125  West  but  are  not  sure  what 
ITD's  you  should  be  receiving  from  9940W  and  9940Y.  To 
determine  these  ITD's,  proceed  as  follows: 


See 

Key. in 

£JL££L£ 

XEQ"PR" 

LAT? 

35 

R/S 

LON? 

125 

R/S 

9940W : 

16019.35 

R/S 

9940Y : 

42584.71 

R/S 

NEXT  OPTION? 

You  should  expect  to  receive  an  ITD  of  16019.35  from  9940W 
and  an  ITD  of  42584.71  from  9940Y. 
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data  in  the  calculator  to  a  known  position  when  the  indicated 
time  delay  (ITD)  is  known  for  each  station  pair. 

Example:  Suppose  you  are  receiving  an  ITD  of  16308  from 
9940W  and  42800  from  9940Y.  These  ITD's  would  tell  you  that 
your  location  is  36d47'55"N  and  121d47 ' 11 "W.  However,  you  know 
that  your  position  is  bench  marked  to  be  at  36d47'36"N  and 
121d46 1 58 "W,  and  you  wish  to  calibrate  your  calculator  so  that 
the  ITD's  of  16308  and  42800  will  give  you  the  latter  fix 
instead  of  the  former.  Proceed  as  follows: 


See 

Key  in 

Press 

XEQ"CA" 

LAT? 

36.4736 

R/S 

LON? 

121.4658 

R/S 

ITD  9940W : 

16308 

R/S 

ITD  9940Y : 

42800 

R/S 

6.  NEXT  OPTION? 


Entering  16308  and  42800  into  the  FI  routine  will  now  give 
you  a  fix  at  36d47'35"N  and  121d46'59"W.  The  small  discrepancy 
between  this  fix  and  the  bench  mark  is  due  to  assumptions  made 
in  the  fixing  algorithm. 

Calibration  is  achieved  by  modifying  the  Master/Slave 
baseline  (BL  in  the  CE  and  ED  routines) .  See  Section  F. 
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The  SW  utility  allows  the  user  to  swap  the  station  data 
stored  in  R16  -  R25  with  the  data  stored  in  R26  -  R35.  Whenever 
the  MI,  XR  or  TR  routine  is  used,  the  SW  routine  is  invoked 
prior  to  the  loading  of  the  data;  the  incoming  data  are  then 
placed  in  R16  -  R25.  The  consequence  is  that  the  first  station 
pair  data  reside  in  R26  -  R35  and  the  second  station  pair  data 
reside  in  R16  -  R25 . 

One  user  application  of  SW  would  be  to  change  the  order  of 
the  station  ID  query  in  the  FI  routine  (this  also  affects  the 
order  of  determination  of  the  solution  and  alternate  solution). 
Another  user  application  would  be  to  output  both  resident 
station  pairs  to  card,  extended  memory  or  tape  cassette  using 
the  CS,  XS  or  TS  routines,  respectively.  To  accomplish  this, 
first  use  the  CS,  XS,  or  TS  routine;  then  XEQ"SW"  (note  that 
the  ID  of  the  station  data  in  R16  -  R25  appears  in  the  display 
instead  of  the  NEXT  OPTION?  prompt) ;  and  finally  use  the  CS,  XS 
or  TS  routine  once  more. 
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11a.  Recording  the  Loran-C  Program  onto  Magnetic  Cards 

(1)  Attach  the  card  reader  to  the  HP-41 CV. 

(2)  Place  the  calculator  in  the  USER  mode. 

(3)  Press  the  PRGM  key  to  place  the  calculator  in  the 
program  mode. 

(4)  Pass  one  side  of  a  blank  magnetic  card  through  the 
card  reader.  Then  follow  the  display  prompts  until 
nine  program  cards  (17  tracks)  have  been  recorded. 

(5)  Press  the  PRGM  key  once  more  to  leave  the  program 
mode. 

(6)  To  record  a  status  card,  XEQ"WSTS".  Then  pass  a  blank 
card  through  the  card  reader  following  the  display 
prompts . 

lib.  Recording,  tiie.  Lgxflfl.-C„£xpgxdffl  on-fco  Magnetic,  Tape 

(1)  Attach  the  tape  cassette  to  the  HP-41CV. 

(2)  Place  the  calculator  in  the  USER  mode. 

(3)  Press  the  alpha  key,  key  in  the  word  LORANC,  press  the 

alpha  key  once  more,  and  then  XEQ"WRTP". 

(4)  To  record  the  program  status,  press  the  alpha  key,  key 

in  the  word  STATUS,  press  the  alpha  key  once  more,  and 

then  XEQ"WRTS" . 
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12a.  Loading  the  Loran-C  Program  from  Magnetic  Cards 


(1)  Attach  the  card  reader  to  the  HP-41CV. 

(2)  Clear  program  memory:  Turn  the  calculator  off,  then, 
while  pressing  the  left  arrow  (erase)  key  down,  turn 
the  calculator  on. 

(3)  Place  the  calculator  in  the  USER  mode. 

(4)  Read  in  the  STATUS  card.  The  status  card  will  set  the 
calculator  to  SIZE  42. 

(5)  Read  in  the  nine  program  cards  (17  tracks). 

12b.  Loading  the  Loran-C  program  from  Tape  Cassette 

(1)  Attach  the  tape  cassette  to  the  HP-41CV. 

(2)  Clear  program  memory.  (See  Step  2  above). 

(3)  Place  the  calculator  in  the  USER  mode. 

(4)  Press  the  alpha  key,  key  in  the  word  STATUS,  press  the 
alpha  key  once  more,  and  then  XEQ"READS".  The  status 
file  will  set  the  calculator  to  SIZE  42. 

(5)  Press  the  alpha  key,  key  in  the  word  LORANC,  press  the 
alpha  key  once  more,  and  then  XEQ"READP". 
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E.  The  Loran-C  Fixing  Algorithm 

The  principles  of  Loran  lines  of  position  (LOP's)  and 
fixing  are  adequately  covered  in  Reference  1  and  will  not  be 
repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 
ITD  =  [Ts  +  p(Ts)]  -  [Tm  +  p(TM)]  +  [TB  +  p(TB)]  +  6  (1) 

where 

ITD  is  the  "indicated  time  difference"  in  microseconds, 

Tw  is  the  distance,  in  microseconds,  from  the  master 

M 

to  the  receiver, 

Tg  is  the  distance,  in  microseconds ,  from  the  slave  to 
the  receiver, 

Td  is  the  distance,  in  microseconds ,  between  the  master 

D 

and  the  slave, 

6  is  the  assigned  station  pair  coding  delay,  in  micro¬ 
seconds,  and 

p (T)  is  the  secondary  phase  correction,  in  microseconds, 
for  a  surface  seawater  path  of  length  T  . 

The  quantity 


At  =  TB  +  p (Tg)  +  6 

is  a  constant  for  each  master/slave  pair.  The  quantity  Tg  is 
computed  from  the  positions  of  the  master  and  slave  using  the 
reverse  solution  algorithm  (Section  H)  at  the  time  of  manual 
data  input  (Routine  MI) . 
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The  following  World  Geodetic  System  1972  (WGS  72)  values 
have  been  adopted  for  Loran-C  navigation  [Ref.  4]: 

Vq  =  299792458  meters/second  is  the  velocity  of  light 
in  free  space, 

n  =  1.000338  is  the  index  of  refraction  of  the 

surface  of  the  earth  for  standard  atmosphere 
and  100  kHz  electromagnetic  waves, 
ag  =  6378135.000  meters  is  the  equatorial  radius  of 
the  earth 

and  f  =  1/298.26  is  the  flattening  factor  (l-b/ae,  where 

b  is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase  cor¬ 
rection  p (T)  are  contained  in  Reference  4,  but  for  use  in  the 
HP-41CV,  the  form 
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On  the  surface  of  a  sphere,  a  hyperbolic  line  of  position 
(LOP)  can  be  represented  by  the  equation  [Ref.  1,  page  175] 


tan  r 


cos  2a  -  cos  2c _ 

sin  2c  cos  w  +  c,  sin  2a 


(3) 


where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c  is  the  spherical  arc  joining  the 
foci,  2a  is  a  constant  for  any  one  hyperbola,  and  r  and  to 
are  the  spherical  coordinates  of  a  point  on  the  hyperbola.  If 
the  base  line  of  the  coordinate  system  is  the  arc  joining  the 
foci  the  to  is  the  spherical  polar  angle  from  the  baseline  to 
a  point  P  on  the  spherical  hyperbola  and  r  is  spherical 
polar  distance  (or  arc)  from  the  prime  focus  to  P  .  Using  the 
Loran  system  we  take  £  =  +1  if  the  prime  focus  is  at  a  master 
station  and  £  =  -1  if  the  prime  focus  is  at  a  slave  station. 

If  we  let  v  =  vQ/n  be  the  velocity  of  100kHz  electro¬ 
magnetic  radiation  at  the  earth's  surface  then,  for  a  spherical 
earth,  we  can  relate  the  parameters  in  Equations  1  and  3  as 
follows : 


2c  =  vTB/ag  , 

and 

2a  =  v(Ts  -  TM)/ae  . 

Using  the  spherical  approximation  for  now,  we  see  that  2c  is 
known  for  any  Loran  pair.  The  "indicated  time  delay"  ITD  is 
measured  by  the  receiver  at  point  P  ,  and  to  determine  a 
hyperbolic  line  of  position  we  must  determine  2a,  but  T„  -  T„ 

B  M 

cannot  be  computed  from  Equations  1  and  2.  If  Sq  were  zero 
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in  Equation  2,  then  it  would  be  possible  to  determine  Tg  -  TM 
uniquely.  As  an  approximation  we  use  the  following  parameters 
in  Equation  2: 


a^  =  -0.321  , 

and  a 2  =  0.000635  . 

These  values  have  been  obtained  by  setting  =  0  and  determin¬ 
ing  a^  and  a2  by  linear  regression  of  the  T  >  537  values  over 
the  interval  of  1000  <  T  <  8000.  This  approximation  is  quite 
good  (within  0.03  ys)  for  distances  up  to  10,000  microseconds 
where  small  changes  in  the  LOP's  can  cause  large  position  errors. 
At  short  distances  the  error  increases  from  0.05  ys  at  1000  ys 
to  0.58  ys  at  10  ys;  although  these  errors  are  large  for  small 
distances,  the  LOP's  are  not  as  sensitive  to  these  changes  as 
they  would  be  at  large  distances.  When  this  approximation  is 
substituted  into  Equation  1,  we  obtain 

[Tg  +  +  a2Tg]  -  [Tm  +  ax  +  a2TM]  =  ITD  -  At  , 

or 

Ts  -  TM  =  (ITD  -  A t)  /  ( 1  +  a2)  (4) 

and  hence  2a  =  v(Tc-TM)/a  is  determined  for  use  in  the 
spherical  approximation. 

Consider  a  Loran-C  triplet  with  the  master  stations  co¬ 
located.  Let  and  C2  denote  the  azimuth  angles  of  slave 

1  (S^)  and  slave  2  (S2> ,  respectively,  measured  from  North 
toward  the  East  from  the  master  stations  (M)  (see  Figure  3) . 
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Further,  let  a  and  r  be  the  azimuth  and  spherical  polar 
arc  (distance)  of  the  receiver  (R)  from  M  .  For  this  geometry, 
Equation  3  can  be  written  in  the  form 


where 


tan  r .  = 


B . 
1 


i  cos  (a  -  £j_)  +  A^  ' 


A.  =  £ .  sin  2a.  , 

l  i  l 


(5) 


and 


B.  =  cos  2a.  -  cos  2c.  , 
i  i  l 


C.  =  sin  2c. 

l  l 


for  the  l  Loran  pair,  i  =  1,2.  Since  r^  =  r2  =  r,  we  can 
eliminate  tan  r  between  the  two  equations.  The  resulting 


equation  can 

be  rewritten  as 

C  cos  a  +  S  sin  a 

=  K  , 

(6) 

where 

C  =  B^C2  cos  ^2  - 

B2C^  cos  £^  , 

S  =  B^C2  sin  £2  “ 

B2C1  sin  £2  , 

and 

K  =  ~  B^A2  . 

If  we  define 

p  >  0  and  y  by 

the  equations 

p  cos  y 

=  c  , 

(7) 

and 

p  sin  y 

=  s  , 

then 

II 

CL 

+  s2  , 

and 

Y  =  qatn 

(S,C)  . 
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Here  the  function  qatn(y,x)  is  the  arctangent  of  y/x  adjusted 
for  the  proper  quadrant  according  to  the  signs  of  x  and  y  . 
A  compact  form  of  this  function  is 


qatn(y,x)  =  tan  1 - q— ^ -  +  nt(x  <  0?) 

x  +  10“y  t ( x  =  0?) 

where  t(z)  =  1  when  z  is  true 

and  t(z)  =  0  when  z  is  false. 

When  convenient  we  will  use  the  notation  qatn(y/x)  inter¬ 
changeably  with  qatn(y,x).  Now  we  can  substitute  Eq.  (7)  into 
Eq.  (6)  and  solve  for 

a  =  y  -  cos  1(K/p)  (8) 

to  obtain  the  azimuth  angle  a  of  the  two  points  of  inter¬ 
section  of  the  spherical  hyperbolic  LOP's.  Finally  we  can 
obtain  a  value  for  r  by  substituting  each  a  into  Eq.  (5) 
for  either  i  =  1  or  i  =  2 .  We  find  that 


r  =  qatn 


B  . 

1 


C  .  cos  (a -E, .)  +  A. 
l  i  l 


for  i  =  1  or  2 , 


(9) 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of  r  . 

The  fixing  algorithm  then  uses  a  and  r  in  the  direct 
solution  algorithm  of  spheroidal  geodesy  (Section  G) . 
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7.  The  Calibration  and  ITD  Prediction  Algorithms 

Calibration  can  be  achieved  when  an  ITD  is  measured  at  a 
known  bench  marked  position.  From  the  bench  marked  position 
and  the  known  master  and  slave  positions,  the  Quantities 
TM  +  P<V  an<^  Ts  +  P<V  can  comPutec^  using  the  reverse 
solution  algorithm  (Section  H)  and  the  accurate  secondary  phase 
correction  formula  (Eq.  2) .  Equation  1  can  then  be  solved  for 
T_  +  p(T_)  to  obtain  a  modified  baseline.  This  modified  base- 
line  is  stored  and  then  used  instead  of  the  true  baseline  in 
subsequent  computations.  The  affect  on  the  accuracy  of  fixes 
using  this  modified  baseline  with  positions  far  removed  from 
the  bench  mark  has  not  studied. 

The  ITD  prediction  algorithm  is  a  direct  application  of 
Equation  1.  A  known  position,  together  with  the  known  master 
and  slave  positions,  is  used  to  compute  the  quantities 
Tm  +  P(TM)  and  Tg  +  p(Tg).  When  these  values,  along  with  a 
computed  or  calibrated  baseline,  TD  +  p(T_)  ,  are  substituted 

into  Equation  1,  a  predicted  ITD  is  obtained. 
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G .  The  Direct  Solution  Algorithm 


This  direct  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8].  Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 

The  qatn  function  is  defined  in  a  previous  section.  East 
longitudes  and  South  latitudes  are  negative.  We  are  given  the 
point  on  the  spheroid,  where  A^  are  the  geodetic 

latitude  and  longitude  (geographic  coordinates) ;  the  forward 
azimuth  a^2  and  distance  S  to  a  second  point  ;  an<^ 

from  these  we  are  to  find  the  geographic  coordinates  <p^f  A 2  and 
the  back  azimuth  a21  .  The  given  quantities  are  i}^,  A^,  a12 
and  S  .  No  assumptions  about  the  relative  location  of  P-j^ 
and  P2  are  required.  The  modified  direct  solution  algorithm 
is : 

0X  =  tan  ^[(l-fjtan  0^],  M  =  -  sin  a^2  cos  0^  , 

C1  =  fM,  C2  =  f ( 1-M2 ) / 4  , 

D  =  (1-C2) (1-C2  -  CXM) ,  P  =  C2[l  +  (1/2)C1M]/D  , 

N  =  cos  0^  cos  a^2  ,  =  qatn(N,  sin  6^)  , 

d  =  S/(aeD)  ,  u  =  2(o^-d)  ,  W  =  1-2P  cos  u  , 

V  =  cos(u+d)  ,  X  =  C2  sin  d  cos  d(2V^-l)  , 

Y  =  2PVW  sin  d,Ao=d+X-Y, 

7  2  1/2 

K  =  [M  +  (N  cos-A a  sin  0^  sin  A a)  ]  , 

0 2  =  tan  ^[(sin  0^  cos  Ao  +  N  sin  Aa)/K]  , 

An  =  qatn(-  sin  Ao  sin  a^2  ,  cos  0^  cos  A  a  -  sin  0^  sin  A  a  cos  a^2)  , 

H  =  C^(l-C2)  Aa  -  C^C2  sin  A a  cos(2a^  -  A a)  , 

a2  =  A1  +  An  -  H  , 
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(*2^  =  qatn[-M,  -  (N  cos  A  a  -  sin  0^  sin  Ao)  ]  , 

<p 2  =  tan-1  [tan  02/(l-f)3 

Details  of  the  modifications  made  to  Thomas'  algorithm  are 
contained  in  Reference  3.  The  algorithm  above  has  been  further 
modified  so  that  Eastern  longitudes,  rather  than  Western 
longitudes,  are  negative. 


H.  The  Reverse  Solution  Algorithm 


This  reverse  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  8-10].  Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The  qatn  function  is  defined  in  a  previous  section.  East 
longitudes  (1)  and  South  latitudes  (4>)  are  negative.  We  are 
given  the  points  P^(<f>^,A^)  ,  P2(4>2,A2)  on  t*ie  sPheroid  and  are 
to  find  the  distance  S  between  the  points  and  the  forward 
and  back  azimuths,  a^2  and  a2^  •  Given  quantities  are 
4>2'  A^,  <p 2  and  A 2  *  No  assumptions  about  the  relative  location 
of  and  P2  are  required.  The  modified  reverse  solution 
algorithm  is: 

6^  =  tan  ^[(l-fltan  (J)^]  ,  i  =  1,2  , 

AA  =  A2  -  ^  ,  A0m  =  (02  -  01)/2  ,  6m  =  ^  +  6 2 ) / 2  , 

H  =  cos2  A0  -  sin2  0  ,  L  =  sin2  A0  +  H  sin2(AA/2)  , 

m  m  m 

d  =  2  sin  1/L  ,  U  =  2  sin2  6  cos2  A©  /(1-L)  , 

m  m 

V  =  2  sin2  A 6  cos2  6/L,X=U+V,Y=U-V, 
mm 

T  =  d/sin  d  ,  D  =  4T2  ,  E  =  2  cos  d  ,  A  =  DE  , 

C  =  T  -  (A-E)/2,  nx  =  X(A+CX)  , 

B  =  2D,  n2  =  Y (B+EY)  ,  n3  =  DXY  , 

<$2d  =  f2(nx  -  n2  +  n3>/64  ,  6jd  =  f(TX-Y)/4  , 

S/ae  =  (T  -  <5^  +  62d)  sin  d  ,  M  =  32T  -  (20T-A)X  -  (B+4)Y  , 

F  =  2Y  -  E ( 4-X)  ,  G  =  f T/2  +  f2M/64  ,  Q  =  -  (FG  tan  AA)/4  , 

AA,!  =  (AA  +  Q)/2  , 
m 
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t,  =  qatn(sin  A0  cos  AA'  ,  cos  0  sin  AA')  , 

1  ^  ra  m  m  m 

t0  =  qatn(-  cos  A0  cos  AA'  ,  sin  0  sin  AA')  , 

2  m  m  m  m 

a12  =  fcl  +  fc2  '  a21  =  fcl  -  fc2 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3.  The  algorithm  above  has  been 
further  modified  so  that  Eastern  longitudes,  rather  than 
Western  longitudes,  are  negative. 
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i.  Program  Accuracy 


The  direct  and  reverse  solution  algorithms  are  equivalent  to 
the  second  order  flattening  algorithms  given  by  Thomas  (Ref.  5); 
the  parameters  of  the  W GS  1972  spheroid  are  used.  The  reverse 
solution  algorithm  reproduces  the  baselines  provided  by  the 
Defense  Mapping  Agency  for  all  40  Loran-C  stations  to  within  0.15 
meters  (the  average  deviation  is  -0.031  meters,  DMA  minus  HP-41CV, 
with  a  standard  deviation  of  0.037  meters)  and  to  within  0.01 
microseconds,  including  the  secondary  phase  correction  for  an  all 
seawater  path.  The  reverse  solution  algorithm  is  also  used  to 
generate  predicted  lTD's;  these  are  presumed  to  be  within  0.01 
microseconds,  also. 

The  fixing  algorithm  uses  the  direct  solution  algorithm  with 
the  azimuth  and  distance  of  the  fix  from  the  vertex  of  the  Loran-C 
triplet  computed  from  Equations  8  and  9  as  inputs.  Equations  8  and 
9  are  based  upon  spherical  geometry  and  include  an  approximation 
to  the  secondary  phase  correction  for  an  all  seawater  path.  The 
largest  source  of  error  is  the  assumption  that  the  azimuth  and 
distance  to  the  fix  are  accurately  represented  by  this  spherical 
approximation.  This  approximation  has  not  been  rigorously  tested, 
however  it  is  possible  to  use  the  reverse  solution  algorithm  to 
predict  the  lTD's  that  will  be  received  at  a  given  position  and 
then  to  enter  these  lTD's  into  the  fixing  algorithm  to  determine 
how  accurately  the  fixing  algorithm  reproduces  the  original 
position.  The  distance  between  the  fix  and  the  original  position 
can  be  determined  using  the  HD  algorithm.  Tables  1,  2  and  3  were 
producted  in  this  manner.  Similiar  tables  with  different  station 
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pairs  are  given  in  References  3,  6  and  7.  It  is  felt  that  the 
results  are  accurate  enough  to  not  warrant  the  inclusion  of  an 
iterative  improvement  routine.  Of  the  samples  in  Tables  1,  2  and 
3,  the  largest  error,  0.22  n.mi.,  is  the  first  entry  in  Table  1. 
From  the  chart  LCNC-2,  it  estimated  that  the  angle  of  intersection 
of  the  two  hyperbolic  lines  of  position  is  about  5  degrees  and  so 
an  error  of  0.22  n.mi.  should  not  be  unexpected. 
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Table  1.  Colocated  Master  Stations 


Position 

Lat  Long 

Predicted  ITD's 

9940W  9940X 

HP-41CV  Fix 

Lat(N)  Long(W) 

Error 

n.mi 

31° 

123° 

16413.28 

27570.93 

30°59’47" 

123600'03" 

0.22 

37° 

126° 

15610.11 

27020.50 

36059'53'' 

126°00' 09" 

0.17 

42° 

129° 

13881.78 

27285.58 

42°00'00" 

129°00'00" 

0 

44° 

132° 

13180.89 

27371.19 

44°00'00" 

132o00'0r' 

0.01 

O 

00 

135° 

12301.25 

27552.06 

48°00'01" 

135°00'03" 

0.03 

50° 

138° 

12068.67 

27584.22 

50°00'01" 

138°00'04" 

0.05 

Table  2.  Colocated  Slave  Stations 


Position 

Lat  Long 

Predicted  ITD's 

9940W  5990Y 

HP- 4 1CV  Fix 

Lat(N)  Long(W) 

Error 

n.mi 

31° 

123° 

16413.28 

27177.18 

3l°00'03" 

122059'58" 

0.06 

37° 

126° 

15610.11 

27403.20 

37°00'01" 

125°59 ' 59" 

0.01 

42° 

129° 

13881.78 

27955.45 

42°00'00" 

129°00'00" 

0 

44° 

132° 

13180.89 

28512.90 

44°00'00" 

132°00'00" 

0 

O 

00 

135° 

12301.25 

29413.61 

48°00'00" 

134°59'59" 

0.01 

Ui 

o 

o 

138° 

12068.67 

29816.84 

50°00'00" 

137°59'59" 

0.02 

Table  3.  Colocated  Master  and  Slave  Stations 


Position 

Lat  Long 

Predicted  ITD's 

5930Y  9960W 

HP-41CV  Fix 

Lat(N)  Long(W) 

Error 

n.mi 

44° 

63° 

29864.46 

11685.15 

44°00'00" 

63°00'00" 

0 

41° 

66° 

30585.61 

12946.91 

41°00'00" 

66°00'00" 

0 

39  0 

69° 

31020.46 

14111.31 

39° 00' 00" 

69°00' 00" 

0 

35  o 

72° 

31064.57 

15139.48 

35° 00' 00" 

72°00’00" 

0 

30° 

75° 

31040.82 

15610.46 

29  °59 1 59" 

75°00*01" 

0.02 

26° 

O 

00 

31106.20 

15858.46 

25059. 57.. 

78600'02" 

0.05 
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Appendix 

A:  Program  Storage  Allocations,  Flag  Usage 

and 

Program  Listing 

Program  Size:  42 

Registers : 

R00  -  R13 : 

Scratch 

storage 

R14 :  Flattening 

2nd  Sta. 

Pair 

Variable 

1st  Sta.  Pair 

R15: 

2a 

R16 : 

ID 

:  R26 

R17 : 

tb  +  p(tb) 

:  R27 

R18 : 

2c 

:  R28 

R19  : 

CD 

:  R29 

R20; 

0  master 

:  R30 

R21: 

A  master 

:  R31 

R22 : 

a  master-slave 

:R32 

R23 : 

0  slave 

:  R33 

R24 : 

A  slave 

:  R34 

R25 : 

a  slave-master 

:  R35 

2a 

:  P36 

R37 : 

0  fix 

R38 : 

A  fix 

R39 : 

a  fix-dest 

R40 : 

0  dest 

R41 : 

A  dest 
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Flag 

Usage : 

FOO : 

1. 

Requery  erroneous  ITD  at  LBL06 

2. 

Echo  two  data  sets  at  LBL  ED 

FOl: 

PR  and  CA  interlock 

F03: 

Set  for  2nd  solution 

F05 : 

Set  if  slave  at  vertex  of  2nd  station 

pair 

F06 : 

Set  if  slave  at  vertex  of  1st  station 

pair 

F07 : 

Set  for  XMEM  functions,  clear  for  tape 

functions 

F08 : 

DN  and  DH  interlock 

F09 : 

Set  if  DMS  conversion  required 

FlO : 

1. 

SW  and  AS  interlock 

2- 

Input  interlock  in  FI  vertex  check 

FI  4: 

PR  and  CA  loop  control 
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Station  U 

o.  of  Pair 

£  Location 

4990 

2 

Central  Pacific 

5930 

2 

East  Coast,  Canada 

5990 

3 

West  Coast,  Canada 

7930 

3 

North  Atlantic 

7960 

2 

Gulf  of  Alaska 

7970 

4 

Norwegian  Sea 

7980 

4 

Southeast  U.S.A. 

7990 

3 

Mediterranean  Sea 

8970 

3 

Great  Lakes 

9940 

3 

West  coast,  U.S.A. 

9960 

4 

Northeast  U.S.A. 

9970 

4 

Northwest  Pacific 

9990 

3 

North  Pacific 
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Appendix  C:  LORAN-C  STATION  DATA 

The  following  list  contains  the  pertinent  parameters  for  each 
Loran-C  station  pair.  This  list  was  compiled  from  the  data  in 
Reference  4.  Each  column  contains  the  following  information : 

1.  The  Loran-C  station  pair  designator. 

2.  The  coding  delay. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  list,  positive  longitudes  are  West,  negative 
longitudes  are  East,  positive  latitudes  are  North  and  negative 
latitudes  are  South.  In  columns  3  through  6  the  latitudes  and 
longitudes  appear  to  be  in  decimal  form,  but  the  actual  format  is 
DDD.MMSSFF  (which  is  compatible  with  the  HP-41CV  D.MS  or  H.MS 
input  mode)  where 

DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds,  and 
FF  designates  hundredths  of  seconds. 
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Appendix  D;  Cfl LOCATED  STATIONS 


5930  MASTER 

— 

9960W 

5930X 

- 

9960X 

5930Y 

- 

7930Z 

5990X 

- 

7960Y 

5990Y 

- 

9940W 

7930W 

- 

7970Y 

7930X 

- 

7970  MASTER 

7960X 

- 

9990Z 

7980Z 

- 

9960Y 

7980  MASTER 

- 

8970W 

8970X 

- 

9960  MASTER 

8970  MASTER 

— 

9960Z 

Slave  stations  are  denoted  with  a  letter 
suffix.  Master  stations  are  so  designated. 
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DISTRIBUTION  LIST 


NO.  OF  COPIES 

Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22314 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Office  of  Research  Administration  1 

Code  01 2 A 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Library,  Code  55  * 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Office  of  Naval  Research  2 

Fleet  Activity  Support  Division 

Code  ONR-230 

800  North  Quincy  Street 

Arlington,  VA  22217 

itevy  Tactical  Support  Activiity  2 

P.G.  Box  1042 

Silver  Springs,  MD  20910 

Attn:  Mary  Ellen  Lannon 

COMPATWINGSPAC  2 

Naval  Air  Station 
Moffett  Field,  CA  94035 
Attn:  Code  51  and  Code  532 

COMPATWINGSLANT  2 

Naval  Air  Station 
Brunswick,  ME  04011 
Attn:  Code  N7 

Prof.  R.H.  Shudde,  Code  55Su  100 

Naval  Postgraduate  School 
Monterey,  CA  93940 
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